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The critical fluctuations in SrTiOa near its 105 K structural phase transition were studied with 
triple axis diffractometry using high energy (> lQOkeV) synchrotron radiation in different SrTiOa 
crystals with different oxygen vacancy concentrations. Due to the presence of oxygen vacancies 
the critical behaviour is changed compared to defect-free systems. In our experiments a smearing 
out of the squared order parameter and a crossover of the critical exponents v and 7 above the 
phase transition temperature is observed, with the crossover temperature strongly depending on the 
concentration of the defects. To understand the experimental findings, e.g. the unusual values for 
the critical exponents found near the critical temperature, the Ginzburg-Landau- Wilson functional 
for structural phase transitions in disordered systems was analyzed using renormalization group 
theory and the replica trick. Considering the effects of defects which locally increase the transition 
temperature leads to a qualitative understanding of the observed behaviour. The crossover behaviour 
of the critical exponents can be modeled and a quantitative analysis of the observed experimental 
data is presented. 
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I. INTRODUCTION 



The influence of defects on the critical behaviour close to a structural phase transition has been studied both 
experimentally and theoretically for a long time. These investigations show that the critical behaviour can change 
dramatically compared to homogeneous systems. According to the Harris' criterium the critical behaviour is modified 
by quenched disorder if the critical exponent v of the correlation length obeys the inequality v > 2/d with d the 
spatial dimensionality of the system G]. As a result the critical exponents near the phase transition temperature 
can change their values 0|| . In addition a smearing out of the specific heat can be observed |i|,|5| . Furthermore the 
central peak, experimentally observed at u> = in the dynamical structure factor above the transition temperature, 
e.g. in SrTiC>3 |£|Q|, can be explained as a consequence of defects (§| 0. In recent years it has been found that 
the critical fluctuations above T c reveal a second length scale in SrTiC>3 as well as in many others systems (ll]-[l]| , a 
feature which has been reviewed by R. Cowley jl4| . This phenomenon was observed in surface-near regions (typically 
ss 100 /im) (l^Jlf| and the origin of the second length scale is considered to be the existence of defects |l7j . 
In the present paper, we concentrate on the bulk behaviour where no sharp component occurs, also the central peak 
is not subject of our investigations, the behaviour of the sharp component will be discussed in a forthcoming paper 
(paper II) |l8|. We have probed the cubic-to-tetragonal phase transition in SrTiC>3 at T c w 100 K with synchrotron 
radiation. The high q-space resolution and the high peak to background ratio obtained with triple-axis diffractometry 
at HASYLAB allows for detailed studies of the critical fluctuations at structural phase transitions . 
The critical scattering above the antiferrodistortive transition of SrTi03 has been studied for samples with different 
oxygen vacancy concentrations. We measured the order parameter and the critical exponents of the correlation 
length and the susceptibility, the oxygen vacancy concentration in the samples was determined independently. The 
experimental findings from the scattering experiments are discussed quantitatively in terms of the influence of defects. 
In the theoretical discussion the static critical behaviour at structural phase transitions is considered with a Ginzburg- 
Landau functional, taking into account the influence of defects through a local increase of the critical temperature 
near the defects. The use of the replica trick makes it possible to investigate the static critical behaviour with the 
help of the renormalization-group theory pl| , p2[ . However, in recent years a possible break down of the symmetry 
of the replica functional (replica-symmetry breaking RSB) is discussed p2|-p5[|. Until now the influence of the free 
energy landscapes or "rare regions", which can occur for defect systems, has not been studied for structural phase 
transitions. Here the possibility of local order parameter clusters must be taken into account. 

The paper is organized as follows: In Sec.|l] we present the experimental details and the results of t he scattering 



experiments, followed by a description of the preparation and characterization of the samples. Sec.|IH| provides the 
theoretical analysis of the experimental data, assuming first the qualitative behaviour and then the order parameter- 
states leading to the minimization of the free energy within the partition function. Finally, in Sec. IV we summarize 
and discuss our results. 



II. EXPERIMENTAL RESULTS 



The scattering experiments have been performed on 5 different samples using 120 and 100 keV synchrotron radiation 
at the undulator beamline PETRA II and the wiggler bcamlinc BW5 at HASYLAB, respectively. All data have been 
collected at 3-axis-diffractometers in Laue geometry at the position of the (511)/2 superlattice reflection. In order to 
get both high incident photon flux and adequate q-space resolution in the scattering plane, annealed Si (311) crystals 
have been used as monochromator and analyzer. The resulting resolution (HWHM) in the scattering plane at BW5 
(PETRA II) was Aq x ^^xlO^ 3 !" 1 (Aq x =l x lO" 3 !" 1 ) in the longitudinal direction and 2xlO- 4 A~ 1 < Aq y < 
2xlO _3 A~ 1 in the transverse direction, depending on the mosaicity of the respective sample. Perpendicular to the 
scattering plane the resolution was of the order of Aq z =lxl0 _1 A _1 . The experimental setup is described in detail 
in Refs. pqjlq ]. The deconvolution of the experimental data was performed as described by Hirota et al. p7| ). 
In order to study the effect of oxygen doping, bulk measurements have been performed on a reduced, an oxidized and 
an as grown Verneuil single crystal. Additionally, a flux-grown sample [281] and an almost perfect crystal, grown by 
means of the float-zone technique pgfl , have been investigated. In Table [j the different sample treatments and growth 
techniques are summarized. The behaviour of the sharp component in these samples will be discussed in paper II 
p8| , in the following we only concentrate on the bulk properties. 

The mosaicity of the Verneuil-grown crystals varied between 30" and 100" which is rather broad compared to the 
almost perfect float-zone grown (l"-7") and flux-grown samples (10"). However, the latter two crystals are slightly 
brownish, perhaps due to iron impurities p0| , in contrast to the transparent as-grown Verneuil samples. Impedance 
measurements have been performed in order to determine the oxygen vacancy concentrations in the different samples 
fcll . Interestingly, the amount of oxygen vacancies is about two orders of magnitude bigger for the crystallographic 
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more perfect samples compared to the as grown and oxidized Verneuil-grown crystals. However, the reduced sample 
shows significant changes: the colour changes from transparent to black and the resistance decreases substantially, i. e. 
the amount of oxygen defects increases about three orders of magnitude. Impedance measurements with microelec- 
trodes of different diameters |^2|,^3| revealed a spatially homogeneous conductivity in the samples I and II. Owing to 
the relation between conductivity and oxygen vacancy concentration [ jsT) it can be concluded that the concentration 
of the oxygen vacancies is distributed homogeneously over the respective samples. Assuming a random distribution 
of these vacancies the mean distance of defects results to d = n^ 1 / 3 with n the defect concentration. This distance 
varies between w 240 A for the lower defect concentrations and » 40 A for the higher concentrations [see Table f|. 
In Fig.[j] the integrated intensity of the superstructure reflection (511)/2, normalized to the extrapolated value at 
T = OK, is plotted versus the reduced temperature r = ^- — 1 for the different Verneuil-grown samples and the 
float-zone grown sample. Below T c the intensity is following a power law with the critical exponent 2(3 = 0.68, which 
is well known for SrTiC>3 Above T c we observe a tail which is dependent on the defect concentrations. The tail 

above the critical temperature can originate from both the critical fluctuations and an additional contribution to the 
squared order parameter caused by defects. This can be interpreted in the way that static order parameter clusters 
exist at temperatures higher than T c . Nevertheless, for the determination of the critical temperatures we followed 
the procedure explained in |6), i.e. at T c — AT half of the intensity at T c + AT has been substracted from the data 
of the integrated intensities of the superstructure reflections, which corrects for the effect of the critical scattering 
contribution to the integrated reflecting power of the superlattice reflection. The values of the critical temperatures, 
which have been found for the different samples, are given in Table H. A shift of the critical temperature to lower 
values with raising defect concentration is visible, which has been observed earlier in Refs. fl35 36 1. 
The inverse correlation length kl{T) and the susceptibility Xh{T) are deduced from the critical scattering at the 
(511)/2 superlattice reflection after deconvolution with the experimental resolution function. The critical scattering 
in the float-zone grown sample has been measured over a much larger temperature interval than those of the other 
samples. Here, a crossover between different values for the critical exponents v and 7 obtained from and Xl(T) 

can be identified [see Fig.|j. These data have been taken in the bulk of the float-zone grown sample but the probed 
volume also included two surfaces (the front and the back side). However, due to the relatively large thickness (12 mm) 
of the sample, the contribution of the sharp component to the signal, originating from a thin layer at the surface, was 
less than 1% close to T c and has been neglected. 

Around t c = 0.11 a change of slope is clearly visible for both kl(T) and Xh(T). This crossover temperature r c is in 
excellent agreement with the value of [see Table [Hj, which is defined as following: Tg is the reduced temperature 
at which the correlation length K£ 1 (r ( i) is equal to the mean distance d of the oxygen vacancies. Above this crossover 
temperature the system can be described within mean-field theory, which leads to exponents v — 0.5 and 7 = 1.0, 
but below the crossover temperature the fitted exponents v = 1.19(4) and 7 = 2.89(4) are much higher than theoret- 
ical values. The non-classical exponents expected for this 3d-Heisenberg system with cubic symmetry, calculated by 
LcGuillou et al. @, are v = 0.705(3) and 7 = 1.386(4). 

In Figure || the inverse correlation length for the four other samples is plotted against the reduced temperature. 
Effective exponents v and 7 have been derived by fits of power laws to each data set without taking into account 
the crossover scenario, because the expected crossover temperature tj, is either too small (« 0.02) or too large (0.20) 
to allow an observation of the crossover within the experimental probed temperature range. The effective exponents 
vary between v « 0.7, 7 « 1.4 (as grown and oxidized Verneuil-grown samples) and v « 1.2, 7 w 2.4 (reduced 
Verneuil-grown sample). It should be emphasized that the scaling relation 7 = (2 — nf)v, with 77 w 0.03 Pqj , still 
holds! Also, the scaling relation 2(3 — dv — 7, which yields (3 = 0.34(6) for the float-zone grown sample, is in excellent 
agreement with the experimental data below T c . 

In the following discussion an attempt is made to understand the microscopic nature of these experimental results. 



III. THEORETICAL DISCUSSION 



In our experimental measurements we found a smearing-out of the squared order parameter near the phase transi- 
tion, which can be explained theoretically by the presence of local order parameter clusters above T c Q]35j|l(J. Defects, 
which locally increase the transition temperature, lead to order parameter clusters above the critical temperature of 
the homogeneous system. For SrTiC>3 we expect a rotation of the oxygen octahedra near the position of oxygen 
vacancies to be responsible for finite order parameter clusters above T c El. These rotations can occur around one 
of the three principal directions. Thus a three dimensional Heisenberg model should be used to model the critical 
behaviour. In Fig.|l| the solid line shows the fit to the experimental data using the power law for the normalized 
integrated intensity versus the reduced temperature I 1 / 2 oc t 34 . Because I 1 / 2 is proportional to the order param- 
eter, the value for the critical exponent (3 is in good agreement with the values obtained for the three dimensional 
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Heisenberg model in the homogeneous case |37|. But, in clear difference to the critical behaviour of homogeneous 
systems is the fact that for higher concentrations of oxygen vacancies a more pronounced tail exists [see Fig.[l] and 
Table ||. Of course the critical fluctuations also contribute to the strength of these tails. However, the contribution 
of the existing order parameter clusters and the dependence of the defect concentration is obvious. Unfortunately, for 
the experimental data investigated here, these two contributions cannot be differentiated adequately. Nevertheless, 
we clearly identified the first clear effect of defects on the critical behaviour. In an additional remark, we want to 
stress that our observation of a finite value of the squared order parameter does not mean that there exists a finite 
bulk value for the order parameter above T c . The local order parameter clusters can be oriented in different directions 
in the order parameter space. If all clusters are oriented in one direction, a finite bulk value for the order parameter 
would result [see Ref. 1011. 

For the float-zone grown crystal a crossover of the critical exponents v and 7 can be clearly observed as another 
effect of the defects [see previous Section], taking place at r » 0.11 [see Fig.||]. This is the temperature range where 
the squared order parameter starts to be noticeable different from zero [see Fig.[lJ. In the general picture of the 
influence of defects on the critical behaviour these findings can be understood as follows. Far away from the critical 
region mean-field exponents are expected, because fluctuations are small Q^^[. These mean- field exponents are found 
indeed for the float-zone grown crystal for r > 0.11 [see Fig.^[. Lowering the temperature (r < 0.11) one has to take 
into account fluctuations and a crossover to non classical critical exponents should occur. Without defects, one would 
expect the critical exponents for the 3D Heisenberg model (e. g. v « 0.705 |B7j). We find near the critical temperature 
for the float-zone grown crystal the value of the exponent v = 1.19 [see Fig.M, which differs substantially from the one 
stated above. Because the correlation length becomes larger near the critical temperature, the critical fluctuations 
can become correlated over distances larger than the mean distance between the defects. In this case the impurities 
can change the critical behaviour and a crossover to a different critical exponent is possible. The Harris criterium 
defines the necessary condition |l]]. If this condition is fulfilled, a dependence of the crossover temperature from the 
defect concentration should be observable. For the crystals probed here, with their different defect concentrations, 
we, indeed, find such a dependence [see Figs.||,|| and Table H. As described above the float-zone grown crystal has 
a crossover temperature equal to Td = 0.11. For all other samples, except the reduced crystal, the characteristic 
temperature Td is very small (< 0.023) and thus the crossover itself is not visible in the experimental data. However, 
the critical exponents, obtained from power law fits to the data-sets in Figure ^, are smaller than v — 1.19 but larger 
than the Heisenberg-value 0.705. The crossover from the Gaussian exponents (e. g. v = 0.5) at high temperatures to 
the "defect-induced" exponents (e. g. v « 1.2) close to the critical temperature should lead to effective exponents v e tt 
with values in the range 0.5 < u e ff < 1.2, as observed. 

For the strongly reduced sample the crossover takes place at very high temperatures (r > 0.2) so that the exponent 
v = 1.2, similar to the value observed in the float-zone grown crystal, is found in the whole temperature range probed 
here. 

This again is in good agreement with the assumption that localized defects are responsible for the change in the 
critical behaviour. In the regime near the critical temperature T c a critical behaviour different from the homogenous 
case, i. e. other critical exponents, should be observable. The critical exponents for systems with quenched disorder 
have been calculated before [Q-^] . Those values can be compared with the critical exponents found in our experiments 
(e.g. in an e-expansion for a 3-component order parameter v — 0.5 + e • 9/64 + ... f» 0.64 with e = 4 — D = 1 in 
three dimensions ||). Unfortunately, we find in our experiments much higher values for v and 7. One important 
question arises immediately: can we understand these findings within the concept of universality or are we facing 
non-universal behaviour? A first hint is that we find the scaling relation 7 = (2 — rj)v still valid for all values of the 
critical exponents for the different crystals near the critical temperature [see Figs.|],[|. As already discussed in the 
last section, the scaling relation for (3 holds, too. 

In order to fully understand the experimental data in terms of the influence of defects on the critical behaviour the 
high values for the critical exponents v and 7 near T c need to be explained. In the next two subsections we try to 
provide a theory which can explain the observed values for the exponents and can model the crossover behaviour. 



A. Defect-functional 



To investigate the questions arising from our experimental findings, the influence of local order parameter clusters, 
occuring above the critical temperature T c , on the critical behaviour is taken into account. At low defect concentrations 
the order parameter clusters are supposed to be well isolated. Therefore, in the mean-field picture multiple local 
minima solutions exist, because the localized clusters can have different orientations in order parameter space fl39| . 
The properties of such ground states were discussed for spin-glasses J|2| , random- ferromagnets |^3| , and recently for 
quantum magnets p4j. For random-ferromagnets Dotsenko et al. p3|,[43[ found that the exponents calculated with the 
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help of the replica-method are unstable with respect to a possible breaking of the replica symmetry. New exponents 
result, different from those obtained by the usual replica- method. Even scenarios for possible instabilities of the theory 
are given. From these studies it is obvious that the states of minimal free energy in the partition function have to be 
considered in detail. In the present paper we sketch the theoretical derivation, for a more detailed discussion see e.g. 

The starting point for the theoretical discussion is a Ginzburg-Landau model for continuous structural phase 
transitions in disordered systems with quenched impurities. Until recently, in the standard approach only fluctuations 
around the homogeneous ground state of the order parameter ip(x) = were considered in the replica-method |p3| . 
We will incorporate fluctuations around the states with the lowest free energy in the partition function. The defect- 
functional used for the investigation of the effects of impurities on the free energy |l0(| reads 



r k k 

^(x);C/] = [ d d x i^((a + t/)(^(x)) 2 + (V^(x)) 2 ) + ^E(^(x)) 2 ] 

J L i=l i=l 



(1) 



Here i = 1 . . . k denotes the component of the order parameter field, a — a'(T — T®) is the harmonic coefficient and 
vanishes at the bulk transition temperature Tj? of the pure system. The coefficient b of the non-linear term has to be 
positive. With the short-range potentials 

u= y £u(x-m D ). (2) 

i D = l 

the effects of randomly distributed defects at the sites Xi D (i = 1 . . . Nr>) are described. They can be considered as 
point like defects. 

Applying the replica-trick [ p2| , the averaging over different defect-distributions can be carried out, only taking into 
account fluctuations around the homogeneous ground state y(x) = 0. The discussion of the resulting fixed-points and 
flow equations can be found in Rcfs. pj|45|]. 

To study the additional effect of possible order parameter clusters we start with the saddle-point approximation 

N u 

- Ap(x) + (a + ^ £/(x - x JD )Mx) + M*) 3 = , (3) 

ir> = l 

which gives the minimum solutions of the free energy. To keep the discussion transparent we choose the number of 
order-parameter components equal to one {k = 1). As far as we consider defects which locally increase the transition 
temperature, local order parameter condensates exist above T c , as demonstrated in earlier studies [ ^9| , ^0[ . In Fig.^J 
a typical order parameter configuration is sketched. Thus in the presence of localized order parameter clusters, a 
more complicated solution of the ground state is found. However, as long as the concentration of the defects and 
the order parameter clusters is small, the value of tp(x) outside these regions decays exponentially. The tail of such 
localized order parameter-clusters can approximately be written as <^/ ofc . ( x ) °c | x _x- e~' x ~ Xi D l/ £ Therefore, 
an approximation of the extremal solution is [^3|j24| l 

JVd 

$W(x;C/) = ^^ ) ^? ofc;0 (x) , (4) 

i=l 

with p — 1 . . . L — 2 N possible solutions because of the Ising symmetry in the order parameter space (k = 1) for a 
given defect distribution U and a\ = ±1. No finite order parameter results for well isolated and non- interacting 
order parameter clusters |10|. The clusters are of course temperature dependent, but we neglect this dependence 
here, assuming that in the temperature range where their contribution is most important the cluster structure does 
not change that much. Now we can write the partition function Z in a way that we consider the global minimum 
solutions $( p ) and fluctuations ip(x) around these solutions 

Z[U] = f D<p(x) jr e -n^H^u] +v (^u] = I Dv{x)e -r[v,u) . z[u] , ( 5) 
J p=i 

where the part of Z with the contributions from the saddle point solutions is written as 
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L 

Z[U] = V e-^^M . e " / (* <P) [x;(7]) V(x) 2 +fc [x;E7] V (x) 3 ] _ ^ 

p=l 

Z and Z depend on the actual defect potential U, because the free energy T does. Thus, we are still investigating 
the effect of frozen disorder. Because the system with defects is no longer translationally invariant, the trace in 
the partition function can no longer be treated in the usual way. The replica-method will be used to deal with the 
problem. 

An additional remark concerning the energy barriers should be made. In our approach the existence of large energy 
barriers between the various saddle-point configurations is assumed. Unlike in glasses these barriers are not infinite. 
But as already pointed out by Narayanan et al. [pij in the thermodynamical limit these barriers tend to infinity. As 
a consequence also the dynamics of the clusters should be very slow. This gives us the possibility to consider just the 
statics and to neglect effects of relaxation and reorientation. 

Performing the disorder average with the help of the replica trick by considering n replicas, the effect of order 
parameter clusters can be incorporated. With the saddle-point solution we find 

^ Di Pl ...D Vn e-TZ = i^W.z n [ Vl ... Vn -U}) u , (7) 

with the replica index a and 

pi ...p„=i 

This is the so called replica partition function for our problem. In the following the additional effect of order 
parameter clusters should be studied, thus the contribution of Z n [U] is our objective. The treatment of this kind of 
replica partition function was studied by Dotsenko et al. |^,^3|. The authors analyzed the problems arising from this 
complicated structure for the solution of the partition function. As a result, they introduced two approximations to 
allow a treatment within the usual renormalization group theory, which will be discussed in the following. 

Firstly, a distribution function for the local order parameter clusters <y3; ofc . ( x ) is presumed Q. With that we 
interpret the solutions arising from the saddle-point approximation as random variables, even though for a given 
distribution of defects these solutions are fixed. However, we want to sum over all possible different distributions, 
hence this treatment is adequate. A Gaussian distribution function f2Sjj 



P[{<piok{*)}] =P[$ (p) (x; U)] = ^ c -^ I ddx y^) 2 -v\ ^) 2 ? 



(.=1 



with <p| ofe (x) as the local order parameter cluster is used. The parameter A depends on the distribution of the 
defect-strength, the temperature, and the coupling constant b. Furthermore, $o( x ) is the global solution <E>( p ) with 
vLfc( x ) = V/ofc o( x )- ^ n thi s approximation A is treated as a constant parameter. The crucial point here is the 
assumption that still is a solution of equation (R). Otherwise new interaction terms would be introduced in the 
functional. 

Now the fluctuations have to be considered. For A small, <I>( p )(x) 2 = $o( x ) 2 + 8$( p \x), the fluctuations (5$( p )(x) 
around the saddle point solutions can be integrated out and it results [E3[ 



L 



P =i 



(x) 2 



(10) 



L _ 

\^ e i/ ^ELi 5,1,(po) ( x ) 2 -U i i ij; =1 »''°)(x)N„(x) ! -2* ( J 



p 1 ...p n =i 



defining A' = (1/A - b/2)- 1 . 

As a second approximation Dotsenko et al. suggested a heuristic ansatz to derive the leading contribution to the 
partition function Z n ||(| . From the solution of the random-energy- models (REM) it is known that a partition function 
like Z n takes its maximum for 1 < xq < n. For the REM it's possible to calculate the parameter xq by extremizing 
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the resulting free energy with respect to xq. Unfortunately, in our problem the integration over fluctuations is done 
in a perturbative approach. This is why we cannot calculate the value of Xq for the problem at hand explicitly. For a 
detailed discussion see Ref. Q . 

Applying the two approximations stated above, a breaking of the symmetry of the replica-functional turns out 
to be significant for the Ginzburg-Landau- Wilson functional discussed here |2^,Q, which is a consequence of the 
fluctuations around the minimal free energy solutions ±$o( x )- The form of the resulting functional was first derived 
by Dotsenko et al. |23|]. Because the resulting form of the functional considered here is equivalent to the one discussed 
in Ref. , we want to refer for all the technical details to Ref. ]23| and Ref. jl3| , especially for the derivation of the 
renormalization group flows and the calculation of the fixed points. To keep our discussion tight, we will focus on the 
interpretation of the results. 

B. Crossover scenarios 

After having shown how the breaking of the replica-symmetry can occur for the defect functional in Eq. (Q), now 
the consequences for the critical exponents are discussed. To analyze the crossover effects we concentrate on the 
effective exponent v e ff of the correlation length as a function of the reduced temperature r. 

From the renormalization group flow eqations, which are calculated to one-loop order p3| , the Gaussian fixed-point 
and the Heisenberg fixed-point occur. They both prove to be unstable. But with the number of order parameter 
components /; in 1 < fc < 4 a third, stable fixed point is possible. For example, the result for the critical exponent v 
is 

1 1 3fc(l - ,tq) 

2 2 16(fc-l)-fc.T (fc + 8) ' 1 ' 

with xq = the usual replica symmetric value of the critical exponent v is obtained ||. Hence the effect of the 
breaking of the replica symmetry (where xq ^ 0) on the critical exponents can be interpreted as a reduction of the 
number of order-parameter components in a subtle manner. 

As mentioned above, the value of xq cannot be derived in the approximation considered here. The experimental 
finding that asymptotically near T c the values of the critical exponent v are very similar, hints to the fact that there 
actually is only one xq for a particular class of defect systems. In Fig.|| a set of possible flows of effective exponent 
v e ff versus the reduced temperature is sketched. Depending on the initial values for the coupling constant b and the 
effective strength of the defects c |^3||l4j], different crossover scenarios can result. Starting from the Gaussian fixed- 
point (large r) the effective exponent v e ff starts with its mean-field value. For smaller r the Heisenberg fixed-point, 
calculated to one-loop order [l/v e ff — 2 — |rp§), is reached, if the initial value for c is small. For larger initial values 
of the defect strength a direct crossover from the mean-field fixed-point to the "defect-induced" fixed-point is found. 
The width of the temperature where the crossover occurs, also varies with the initial values of c and b. Assuming 
xq = 0.9622, we can obtain the experimentally found critical exponent v e ff = 1.19 for r — > 0, that is in the critical 
regime. The temperature scale where the crossover takes place is dependend on the initial values of the coupling 
constants. By matching the temperature scale pfj], the crossover can be modeled, too. For the float-zone grown 
sample the resulting fit of the experimental data is shown in Fig.|(| 

Despite the remarkable agreement between theory and experimental findings one should be aware of the fact that 
this is only a one loop result. With a higher order in a loop expansion a similar fit might well require a different value 
for xq. Through the two approximations, defined in the last subsection, a lot of physical effects which are very hard 
to model are condensed into the parameter xq. The next goal of the theory should be to determine the value of :ro 
and to answer the question if this parameter depends on the defect strength or the defect concentration. Dotsenko 
et al. expressed their hope that x$ may be determined in higher order in e, in which case universal exponents for a 
paticular class of defect sysems would result. 

IV. DISCUSSION 

In summary, we have investigated the effects of defects on the critical behaviour at structural phase transitions. 
Our experiments probed the 105 K phase transition of SrTi03 and the influence of oxygen vacancies for different 
SrTiO-3 crystals. 

We found in our experiments that the squared order parameter is non-vanishing and rounded above the critical 
temperature T c , representing the transtion temperature for a perfect crystal. The rounded squared order parameter 
increases above this temperature with increasing concentration of oxygen vacancies in the crystals. 
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These effects were interpreted as a result of local order parameter clusters occuring above T c (see Ref. [[L0| ) . 

Furthermore, near T c unusual high values for the critical exponents v und 7 were found in our experiments. For 
the float zone grown crystal we were also able to study a crossover of the critical exponents v and 7 in great detail. 
In addition, we identified a dependence of the crossover-temperature on the defect concentration of the respective 
samples. 

The crossover scenario was analyzed in the framework of the renormalization-group theory (see also Refs. fi^ , [46|] ). 
Two main effects could be studied here. Firstly, the correlation length of the critical fluctuations has to exceed the 
length scale of the mean distance between oxygen vacancies in order to change the values of the critical exponents 
of the homogeneous system. Secondly, the unusual high values for the critical exponents found experimentally below 
the crossover temperature, which indicate the crossover from classical to non-classical critical behaviour, could be 
explained by taking into account the fluctuations around the states of minimal free energy and the resulting breaking 
of the replica symmetry. 

Because of the occurence of the paramter xo within the theory, which cannot be determined explicitly, some 
questions concerning the universal aspects of the results are still open. Two scenarios are possible: (i) the parameter 
xq is universal for certain classes of defect systems; (ii) the parameter xq depends on microscopic properties of the 
crystals and is not universal. We stated some hints that for the class of defect systems studied here universal behaviour 
occurs. 

It is known from defect models Jitj| that in three dimensions there is a minimum strength of defects. Below this 
border no order parameter clusters form and the analyzed crossover should not occur. 

It would be interesting to study similar systems in the presence of different kinds of defects or to study systems 
where the universality class can be changed, e.g. under pressure. Also further theoretical investigations are necessary 
to get a step beyond the approximations considered in this and in other papers that take into account the breaking 
of the replica-symmetry. 
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Sample ^ 


Growth technique 


sample preparation 


sample colour 


i 


Floatzone grown 




Brownish, transp. 


ii 


Flux grown 


etched (89% H3PO4, lh, 160° C) 


Brownish, transp. 




Verneuil grown 






in 


Oxidized 


(48h, 2 -atm., lbar, 1000° C) 


Rose, transp. 


IV 


As grown 




transparent 


V 


Reduced 


(5h, H 2 -atm., lbar, 1250°C) 


Black 



TABLE I. Specifications of the samples investigated. Specimens grown by three different techniques (bold type) were 
available. 





T c 


n [cm 3 ] 


d[A] 


1 [io-^A- 1 ] 




Floatzone grown 


98.8(1) 


6.1(2) xl0 IS 


55(1) 


18.2(2) 


~ 0.115 


Flux grown 


102.6(2) 


2.8(2) xl0 IS 


71(2) 


14.2(3) 


~ 0.035 


Verneuil grown 












Oxidized 


105.7(1) 


7.4(2) xlO 16 


238(2) 


4.2(1) 


~ 0.014 


As grown 


105.8(1) 


7.6(2) xlO 16 


236(2) 


4.2(1) 


~ 0.025 


Reduced 


101.0(1) 


1.7(1) xlO 19 


39(1) 


25.7(5) 


- 0.20 



TABLE II. Critical temperatures T c and defect concentrations n of the different samples, d = n- 1 ' 3 is the mean distance of 
defects, Td is taken from the measurements of the inverse correlation length («(rd) = d _1 ). 
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Vern. reduced 
Float-zone 
Vern. as grown 
Vern. oxidized 
P = 0.34 




-0.05 0.00 0.05 

reduced temperature x 



0.10 



FIG. 1. Integrated intensities of the superlattice reflection (511)/2. In order to compare the tails above T c the intensities have 
been normalized such that they coincide below the critical temperature. The critical exponent f3 has been fixed to (3 = 0.34. 
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0.01 0.1 1 0.01 0.1 1 

red. temperature t red. temperature t 



FIG. 2. Inverse correlation length and susceptibility of the broad component in the bulk of the Float-zone grown sample. 
Around r c = 0.11 a crossover in the critical behaviour is clearly visible. The inset shows a transverse scan and the best fit to 
the data at T = 100.8 K. 
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oxidized, v=0.73, y=1 .49 
as grown, v=0.79, y=1 .58 
• str. reduced, v=1. 18, y=2.45 



0.010 0.100 
reduced temperature x 



1.000 



FIG. 3. Inverse correlation length for different samples. The effective exponents v and 7 have been derived from power law 
fits to the inverse correlation length and the susceptibility. 
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FIG. 4. A typical order parameter configuration above the critical temperature is sketched. Thus, in the presence of localized 
order parameter clusters a complicated solution of the ground state results. 
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FIG. 5. The effective exponent v is plotted against the temperature-scale for different initial values of the coupling constants. 
Xq is choosen as 0.9622. 
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FIG. 6. Inverse correlation length of the broad component in the bulk of the float-zone grown sample. The theoretical curve 
models the crossover. 
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